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We investigate theoretically the electron spin relaxation in single-electron nanowire-based semicon- 
ductor double quantum dots induced by confined phonons and find that the electron spin relaxation 
rate can be efficiently manipulated by external electric fields in such system. An anti-crossing, due 
to the coaction of the electric field, the magnetic field and the spin-orbit coupling, exists between 
the lowest two excited states. Both energies and spins of the electron states can be efficiently tuned 
by the electric field around the anti-crossing point. Multiple sharp peaks exist in the electric-field 
dependence of the spin relaxation rate induced by the confined phonons, which can be ascribed to 
the large density of states of the confined phonons at the van Hove singularities. This feature sug- 
gests that the nanowire-based double quantum dots can be used as electric tunable on-and-off spin 
switches, which are more sensitive and fiexible than the ones based on quantum-well based double 
quantum dots. The temperature dependence of the spin relaxation rate at the anti-crossing point 
are calculated and a smooth peak, indicating the importance of the contribution of the off-diagonal 
elements of the density matrix to the spin relaxation, is observed. 

PACS numbers: 72.25.Rb, 73. 21. La, 63.22.-m, 63.20.kd 



I. INTRODUCTION 

Manipulation of spin relaxation in semiconductor 
quantum dots (QDs) is of great interest due to its poten- 
tial applications in quantum information processing and 
spintronic devicesj^^— Various mechanisms can lead to the 
spin relaxation, such as the spin-orbit coupling (SOC) 
due to the space inversion asymmetry together with the 
electron-phonon scattering;^"— the coaction of the hyper- 
fine interaction and the electron-phonon scattering^i2r2i 
the g-factor fluctuations and the direct spin-phonon scat- 
tering due to the phonon-induced strain<^ For the widely- 
investigated III-V semiconductor QDs, the SOC together 
with the electron-phonon scattering dominates the spin 
relaxation in large magnetic fields j^^—ii^ii^ It has been 
shown that the spin relaxation can be manipulated not 
only in single QDs by external magnetic field j^^'^^ but 
also in double quantum dots (DQDs) by external electric 
fieldJ^iii The latter is more interesting since the electric 
field is more easily accessible and controllable in gen- 
uine devices. The electric coherent spin manipulation in 
DQDs has been demonstrated experimentally by electric 
dipole spin resonance with weak electric field.-- "^^ It was 
also suggested that the electron spin in DQDs can be 
controlled by strong pumped electric field.-^- 

The above mentioned QDs are either self- assembled 
ones or fabricated by confining electrons in quantum 
wells, where bulk phonons in substrates play an impor- 
tant rolcj^ii^ In recent years, a new type of QDs based 
on self-assembled nanowires has been fabricatedi^"— 
The nanowires are perpendicular to the substrates, mak- 
ing the bulk phonons in the substrates less important 
than the quasi-one-dimensional confined phonons in the 
nanowires i^Si^ It has been shown that the confined 



phonons can lead to novel properties in optical absorp- 
tion in nanowire-based single QDs and transport through 
nanowire-based DQDs<^"— Specifically, the spin relax- 
ation in nanowire-based single QDs can be efficiently 
manipulated by external magnetic field due to the large 
phonon density of states (DOS) near the van Hove singu- 
larity of the confined phonons^^i However, the manipula- 
tion of spin relaxation by electric field in these nanowire- 
based QDs is not reported yet. 

In this paper, we propose a scheme of manipulation of 
the spin relaxation in nanowire-based elongate DQDs by 
external electric field. In our system, the electric field 
can efficiently tune the energy spectrum of the DQDs 
where an anti-crossing between electron states with op- 
posite majority spin exists. As the energy splittings be- 
tween these states match certain van Hove singularity 
of the confined phonons, multiple sharp peaks occur in 
the electric field dependence of the spin relaxation rate 
(SRR). This feature offers an efficient scheme for the ma- 
nipulation of the spin relaxation and hence on-and-off 
switches. Comparing to the electric spin relaxation ma- 
nipulation in quantum-well-based DQDs^^ the SRR in 
our system is much more sensitive to the electric field 
and can be manipulated with more flexibility. Moreover, 
due to the strong spin mixing in the vicinity of the anti- 
crossing point, the off-diagonal density matrix elements 
have a large contribution to the SRR. Such contribu- 
tion can have a pronounced impact on the SRR, which 
can be seen in the temperature dependence of the SRR. 
It should be noted that the strong spin mixing makes 
it difficult to define the majority spin for each state in 
the vicinity of the anti-crossing point. This makes the 
widely-used Fermi-golden-rule approacb^i^rJi for spin 
relaxation calculation become obscure in this regime. So 
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in the calculation, we apply the equation-of-motion ap- 
proach developed by Jiang et al.M^ 

We organize the paper as follows. In Sec. II, we in- 
troduce the model for our system and specify the Hamil- 
tonian. The equation-of-motion approach is also briefly 
described in this section. In Sec. Ill, we first present the 
electric field dependence of the energy levels and the cor- 
responding spin states of the DQDs. Then we show the 
electric field dependence of the SRR and demonstrate the 
electric manipulation of spin relaxation. Finally, we dis- 
cuss the impact of the off-diagonal elements of the density 
matrix on the temperature dependence of the SRR. The 
validation of the Fermi-golden-rule is also addressed in 
this section. We summarize in Sec. IV. 



II. MODEL AND FORMALISM 

We consider a single-electron elongate DQDs embed- 
ded in In As cylindrical nanowire with radius R. An ex- 
ternal magnetic field B is applied along the wire (the 
z-axis). The total Hamiltonian can be expressed as 

Htot — Hc + Hph + Hep, (1) 

where He, i?ph and H^p are the electron, phonon and 
electron-phonon interaction Hamiltonian, respectively. 
The electron Hamiltonian has the form — Hq + Hgo- 
Here Hq is electron Hamiltonian without the SOC: 

i^o = ^ + V,ir) + V,{z) + HB, (2) 
2m* 

with m* denoting the electron effective mass. Vc{r) — 
^m*ujQr'^ is the QD confinement in radial direction. Thus 
the dot has an effective diameter do = ^Jhn /m*ijjQ. The 
confinement Vz{z) in the z-axis is given by 

{eEz + \e{a + 2d)E, \a < \z\ < \a + d, 

eEz + \e{a + 2d)E + Vo, \z\ < \a, 
oo, otherwise, 

(3) 

with Vq representing the barrier height between the two 
dots, a standing for the interdot distance and d denot- 
ing the dot length. A schematic of Vz{z) is plotted in 
Fig. [irb)(IV). In our model, the radial confinement is 
much stronger than the longitudinal one, so we consider 
only the lowest electron subband in the radial direction. 
Hb = ^gi^sB ■ cr is the Zeeman term with g, and 
(T representing the electron g-factor, Bohr magneton and 
Pauli matrix, respectively. The SOC Hamiltonian takes 
the form 

Hso = ■^^yP^' 

where 7 is the Rashba strength .-S> 2^ 

The confined phonon modes are calculated within 
isotropic elastic continuum model by assuming free- 
surface boundary conditions 1^"— The Hamiltonian is 



given by Hph = J2nX ^q\aq\aq\, where w^a is the eigen- 
frequency of the Atn phonon mode with axial wave vector 
q, which is solely determined by the longitudinal or trans- 
verse sound velocity or vt, respectively. Note that the 
calculated confined phonon modes can be classified into 
axial and radial modes according to the corresponding 
displacement field at small g. 30-42 

Both the deformation potential and piezoelectric cou- 
plings are considered in our calculation. Assuming the 
displacement field of the confined phonons is u{r), the 
deformation potential coupling is given by the diver- 
gence of u{r): H^p{r) = — SV • u{r), where S is the 
deformation coupling strength. The piezoelectric cou- 
pling is determined via i/^(r) = ^ f dr^ j'^ | , where 
K denotes the relative dielectric constant, e is the elec- 
tron charge and stands for the polarization induced 
by u{r)^ For InAs nanowires with wurtzite structure, 
pp^(r) = eir,[dzUr{r) + drUz{r)]er + [631(9^ + l)ur{r) + 
Si3dzUz{r)]ez, where 615, 631 and 633 are the piezoelectric 
constants4^i^ By substituting u{r) into the expression 
of and H^p,^'^ one can express the Hamiltonian of 
the electron-phonon interaction in the form 

^op= E E^^^^aX^aKa + (5) 

where = KE? / {2tt pmi^qxR^) for the deformation 

potential coupling and = Sfnre'^el^/ {k'^ R'^cjqx) for 

the piezoelectric coupling, with quantity p,„ being the 
mass density. The factor Xqxi'"') depends on the confined 
phonon eigenmode. Note that since in the calculation, 
we consider only the lowest radial subband of the elec- 
tron, only the dilatation phonon modes couple to the 
electron.— 1^ 

Given the electron Hamiltonian H^, the eigenstates 
{|^)} and corresponding eigenenergies {eg} can be ob- 
tained by exact-diagonalization method.^ The spin of 
each state can be calculated as {(^\S\£) with S — ^cr. 
We apply the equation-of-motion approach to calculate 
the spin relaxation time. Following Ref. [l^ we use the 
density matrix p = p'^ €5 p^^ to describe the system. Its 
time evolution in the interaction picture reads 

dtp\t) = -^[v\t),pm 

-^iyr[V\t),[V\r),p\r)]]. (6) 

By applying the Born-Markovian approximation and car- 
rying out the integral over the constant energy surface of 
the confined phonons, the equation-of-motion of the den- 
sity matrix p*^ = trph(p) for electron can be obtained 
from the above equation within the basis of the electron 
eigenstates (in Schrodinger picture) ^"^ 

dtpl^n - -i ^"^ E E{^mfcfe;(e;fc)pFn 

r]=D,P kl 

-Wl^lu{^ln)pll+^.c], (7) 
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I A,f»7 \2^.Vmk nln 



X {[n{ujgx) + l]0{eni) 
+ n{ujqx)0{ein) 



(8) 



in which Xg™" = (™Ix^'aI"') ^^"^ ^(e/n) denotes the step 
function with £;„ = £; — £„. Note that the transition rate 
r^^^^ from state \^^) to |£2) can be expressed as T^^^^ ~ 

By applying a proper initial state, one can use the 
equation-of-motion of the density matrix to calculate the 
time evolution of {Sz) = tr[Szp''(t)], which is the ex- 
pectation value of the ^-component of the electron spin. 
We assume that the initial state is prepared by electri- 
cal spin injection^i^ or optical pumping of the electron 
spin,— which can polarize the electron spin along the 
external magnetic field. Thus density matrix elements in 
the basis of {|^)} for the initial state can be expressed as 
P«'(0) = where {\Cj} are the spin-down 

eigenstates ((CI<S'z|0 = ^1/2) of with eigenenergies 
{e{\. = Cexp[— £^/(fcBr)] is the Maxwell-Boltzmann 
distribution with C denoting the normalization param- 
eter. Note that due to the SOC, p'^{0) can have off- 
diagonal elements. 

Given the time evolution of {Sz), we define the effective 
SRR as 



= (lna)/Tc 



(9) 



with Ta being the time that the envelope of (Sz) de- 
cays to the 1/a from its initial value to its equilibrium 
value. Note that comparing to the SRR used in previous 
works , ^'^1^° we introduce an additional factor Ina in or- 
der to guarantee that if (Sz) can be well-described by a 
single exponential decay, f ~^ is independent on a. This 
is helpful when we compare the effective SRR with 
different a. It should be emphasized that if (Sz) does not 
follow a single exponential decay, the effective SRR 
is a-dependent. One should properly choose a to give an 
overall description of the spin relaxation in the regime of 
investigation. 



III. NUMERICAL RESULTS 

In the computation, we consider a DQD with length 
d = 50 nm for each dot, interdot distance a = 10 nm, 
barrier height Vq = 0.3 eV and effective dot diame- 
ter do = 12 nm. The Rashba strength 7 is set to 
be 2.4 X 10^^^ eV-m, corresponding to the spin-orbit 
length Ago ~ /{'m*^) — 133 nm which is close to the 
value Aso ^ 127 nm reported in Ref. [2^. The radius 
of nanowire R is set to be 30 nm. The InAs nanowires 
have wurtzite crystal structure with orientation along the 



[0001] directioUi ^^i^^ The piezoelectric constants for the 
wurtzite InAs can be expressed as 615 = 631 = —633/2 = 
— ei4/\/3 with ei4 being the piezoelectric constant for 
bulk InAs.*^ Note that although both the deformation 
potential and piezoelectric couplings are considered in 
our computation, we find that the piezoelectric coupling 
is dominant in our system. The other parameters used 
in our calculation are listed in Table I. 

TABLE I. Material parameters used in the computation 
(from Ref. 52 unless otherwise specified). 



Pm 5900 kg/m^ 

2130 m/s 
vl 4410 m/s 

ei4 3.5 X 10*^ V/r 



m 

K 



5.8 eV 
0.0239 mo 
15.15 

-9.0^ 



Ref. 25. 



A. Energy levels and spin configurations 

Before we study the behavior of the effective SRR, it is 
helpful to first discuss the electric field dependence of en- 
ergies and spins of electron states in the DQD. We plot 
the lowest three energy levels of the DQD as function 
of electric field in Fig. [Ija) with external magnetic field 
B = 1 T. An anti-crossing between the first excited state 
1 2) and the second excited state |3) can be clearly identi- 
fied in the energy levels which is induced by the coaction 
of the electric field, the magnetic field and the SOC. The 
position of the anti-crossing point depends on the mag- 
netic field, which is shown in Fig. [ija) by dots at some 
selected magnetic fields. The arrows in Fig. [ija) illus- 
trate the spin of each state along the z-axis defined as 
{£\Sz\() (The other spin components are zero). With the 
increase of the electric field, the spins of the lowest three 
levels (|1), |2) and |3)) change from f, t and I (before the 
anti-crossing point) to t, i and f (after the anti-crossing 
point), respectively. In the vicinity of the anti-crossing 
point, due to the spin mixing induced by the SOC, the 
spins of the states |2) and |3) decrease and tend to zero 
at the anti-crossing point, while the spin of |1) is not 
sensitive to the electric field. 

The electron spin states can be better illustrated by 
the spin configuration of each state, e.g., the distribu- 
tion of each spin state as function of z. We plot the spin 
configurations of the three states in Fig. [Ijb) at three 
typical electric fields Ei/Eu/Em corresponding to the 
fields before/at /after the anti-crossing point. The com- 
ponent (Sy) of the spin configuration is always zero for 
all the three states so it is not plotted in the figure. One 
finds that the spin configuration of the ground state is 
not sensitive to the electric field and it is mainly dis- 
tributed in the left dot for not too small electric field 
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FIG. 1: (Color online) (a) The lowest three energy levels vs 
the electric field E with the magnetic field B = 1 T. The 
arrows indicate the spins of the states. Longer arrow means 
larger spin. The thin pink dashed curve with dots represents 
the position of the anti-crossing point between the states |2) 
and 1 3) at different magnetic fields, where the corresponding 
magnetic fields of the dots are 0.27 T, 0.46 T, 0.81 T, 1.0 T, 
1.16 T, 1.53 T and 1.9 T, respectively (from left to right). 
The three vertical lines label the typical electric fields Ei = 
0.05 kV/cm, En = 0.156 kV/cm and Em = 0.22 kV/cm 
corresponding to the regime before, at and after the anti- 
crossing point, (b) The spin configurations of the lowest three 
energy levels as function of z for the three typical electric fields 
El, Ell and -Em. Arrow in (b) and curve in (a) with the same 
color and line type stand for the same state. The schematic 
of the potential of the DQD is shown in (b)(IV). 



{E > 0.01 kV/cm). The effect of the electric field on 
the spin configuration is pronounced for the two excited 
states. There are two types of spin configurations for the 
excited states far away from the anti-crossing point: (i) 
the spin configuration distributes mainly in the right dot 
and has a similar structure as the spin configuration of 
the ground state in the left dot; (ii) the spin configura- 
tion distributes mainly in the left dot and has a different 



structure as the ground state. Before the anti-crossing 
point [Fig.[Ii;b)(Ei)], |2)/|3) belongs to type (i)/(ii). In 
contrast, after the anti-crossing point [Fig. [IIb)(Eni)], 
1 2) /1 3) belongs to type (ii) / (i). In the vicinity of the anti- 
crossing point [Fig.[T][b)(En)], both excited states belong 
to a new type of spin configuration, which is evenly dis- 
tributed in the DQD and polarized mainly along the x- 
axis. This suggests that there is a strong spin mixing 
between the two excited states, which makes the expec- 
tation value of the electron spin (5*2) associate with not 
only the diagonal elements but also the off-diagonal ele- 
ments between the two states in the density matrix. 

As both the energies and spins of the two excited states 
can be efficiently manipulated in the vicinity of the anti- 
crossing point, one would expect that they should have 
a pronounced impact on the spin relaxation. So in the 
following discussion, we will concentrate on the lowest 
three levels. 



B. Effective spin relaxation rate 

With the equation-of-motion approach we applied 
here, the effective SRR is obtained from the time evo- 
lution of {Sz), which has different behaviors under dif- 
ferent electric/magnetic fields and temperatures. So we 
start our discussion of the spin relaxation by studying the 
typical time evolution of (Sz)- Then we study the electric 
dependence of the effective SRR and demonstrate how 
the electric spin relaxation manipulation can be achieved 
in our system. We also compare our scheme to the scheme 
based on DQDs in quantum wells. Finally, by studying 
the temperature dependence, we also show the impor- 
tance of the off-diagonal elements of the density matrix 
to the spin relaxation. 

1. Typical behavior of (Sz) 

We show the typical behaviors of time evolution of (Sz) 
in both the short- and long-time regimes at different elec- 
tric fields in Fig.[2Ia)~(d). Generally speaking, they can 
be classified into two types: (i) single-exponential decay 
[Fig. [UJa)] and (ii) double-exponential decay [Fig. [DJb), 
(c) and (d)]. For type (i), the time evolution of (Sz) [dots 
in Fig. mja)] can be well-fitted by a single-exponential 
function fi{t) — aie~*/'^'=i -|- bi [blue dashed curve in 
Fig. mja)] with oi, bi and being fitting parameters. 
For type (ii), the time evolution of (Sz) [dots in Fig.[5](b), 
(c) and (d)] has two stages: a short-time stage where (Sz) 
decays faster and a long-time stage where (Sz) decays 
slower. The relative importance of the two stages can be 
different. For example, the short-time stage is important 
in Fig.[5Jb) since it leads to a large variation of {Sz) from 
—0.5 to 0.3. While in Fig. [2jd), the long-time stage is 
important. Note that usually the time evolution of (Sz) 
in the short-time stage is much faster than that in the 
long-time stage. It is also worth noting that a strong 
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FIG. 2: (Color online) (a) ~ (d) Time evolution of (S^) 
in both the short- (lower scale) and long-time (upper scale) 
regimes with B = 1 T and T = K at the electric field 
E = 0.05, 0.14, 0.146 and 0.16 kV/cm, respectively. Black 
dashed curve with dots: (Sz) from the equation-of-motion 
approach. Blue dashed curve: best-fitting with single expo- 
nential decay. Red solid curve: best-fitting with double ex- 
ponential decay, (e) SRRs as function of electric field with 
B = 1 T and T = K. Red solid/Green dashed curve: ef- 
fective SRR T~^/f^^. Blue dots/Pink triangles: rate for the 
first /second spin relaxation process. The cyan dashed box 
indicates the vicinity of the anti-crossing point. The gray 
vertical line indicates the position of the anti-crossing point. 



oscillation can exist in the short-time stage, as shown in 
Fig.m^c). The evolution of (Sz) can be fitted by a double 
exponential function f2{t) — a\e~^l'^'^^ -t- a2e~*/'^'=^2 -|- h 
[orange-red solid curves in Fig. [51Jb)~(d)] with ai, 02, 6, 
tei and t^.^ being fitting parameters. 

The behavior of {S^ suggests that there are two main 
spin relaxation processes in the system: the first process 
with the characteristic rate t^'^ , which dominates in the 
short-time stage and the second process with the charac- 
teristic rate , which dominates in the long-time stage. 
If the time evolution of (Sz) belongs to type (i), only one 
process dominates (corresponding to t-^^), which solely 
determines the effective SRRs. If the time evolution of 
(Sz) belongs to type (ii), both processes can contribute 
to the effective SRRs. The effective SRR is given 
by the average over the two rates r^T^ and t-^^ with the 



weight determined by the parameter a. For small a, the 
effective SRR mainly describes the short time behav- 
ior and hence has a large weig ht ofT^y While for large 
a, makes a large contribution to the f^^ . 

To give a comprehensive description of both the short- 
and long-time behaviors of the spin relaxation, we choose 
the effective SRR with two typical values of a: a — 
1.5 which better describes the short-time behavior and 
a — e which better describes the long-time behavior. 



2. Electric field dependence 

The effective SRRs fj~g and as function of electric 
field at magnetic field i? = 1 T and temperature T = K 
are plotted in Fig. [UJe). One can see that both effec- 
tive SRRs have similar electric field dependence with two 
main features: one is the large drop at the anti-crossing 
point, the other is the multiple sharp peaks indicated by 
the numbered arrows. 

To understand these features, it is helpful to first 
identify the dominant spin relaxation process at dif- 
ferent electric fields. This can be done by comparing 
the effective SRRs f^l and to the two rates r^T^ 
and ^ in Fig. [H^e). Far from the anti-crossing point 
{E < 0.102 kV/cm and E > 0.19 kV/cm), only one 
spin relaxation process dominates. Both effective SRRs 
are determined by the corresponding rate t-^^ ■ For the 
electric field E in the vicinity of the anti-crossing point 
[E e (0.102,0.19) kV/cm], the second spin relaxation 
process becomes important and can dominate the effec- 
tive SRRs after the anti-crossing point. The difference of 
the two effective SRRs fj"^ and f"-'^ at the anti-crossing 
point is due to the fact that the rate t-^^ has a smaller 
contribution to than to f^^- 

Having identified the dominant spin relaxation process, 
one finds that the first feature, e.g., the large drop at the 
anti-crossing point, is due to the relative importance of 
the two spin relaxation processes: the second spin relax- 
ation process becomes more and more important as the 
electric field E passes the anti-crossing point. The second 
feature, e.g., the multiple sharp peaks, is solely induced 
by the spin relaxation process corresponding to TjT ^ . 

The relative importance of the two processes can be 
understood by studying the relevant energy levels. This 
can be done by comparing the transition rates between 
the lowest three energy levels to the rates t^^ and r^T^ in 
Fig. El^a) . The first spin relaxation process correspond- 
ing to TjT^ is determined by different transitions at dif- 
ferent electric fields. For E < 0.19 kV/cm, this pro- 
cess is mainly determined by the transition |3) — > |1) or 
1 3) — ?> |2) . The transition with larger transition rate dom- 
inates this process. For > 0.19 kV/cm, this process is 
determined by the transition |2) — > |1). The second spin 
relaxation process corresponding to t-^^ , which only ex- 
ists in the regime E G (0.102, 0.19) kV/cm, is determined 
by the transition |2) — ^ |1). So the dominant spin relax- 
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FIG. 3: (Color online) (a) Transition rates (yellow solid 
curve), (purple dashed curve) and Ffi (green dotted 

curve) due to the piezoelectric coupling only, with B = 
1 T and T = K. Blue dots/Pink triangles: rate for the 
first /second spin relaxation process. Inset: the enlargement 
of the peak 4. The cyan dashed box indicates the vicinity of 
the anti-crossing point and the gray vertical line indicates the 
position of the anti-crossing point, which are the same with 
those in Fig. [2jb). (b) Spectrum of the confined phonons 
in the nanowire. The orange solid/blue dashed curves repre- 
sent the axial/radial phonon modes, (c) DOS of the confined 
phonons. The arrows with numbers in (b) and (c) indicate 
the van Hove singularities which are responsible for the sharp 
peaks labelled with the same numbers in (a). 



ation process changes from |3) — > |2) to |2) |1) as 
the electric field E passes the anti-crossing point, which 
agrees with the spins of the states shown in Fig. [Ija). 
This explains the large drop at the anti-crossing point. 
Note that in the calculation of the transition rate, we con- 
sider only the piezoelectric coupling since it is dominant 
in our system.— 

Now let us turn to the second feature of the effective 
SRRs, e.g., the multiple sharp peaks, which are deter- 
mined by the first spin relaxation process corresponding 
to TjT^. These sharp peaks occur as the energy split- 
ting between the two states of the spin relaxation pro- 
cess matches the van Hove singularities of the confined 
phonons, e.g., the denominator of the transition rate 
r^^£^ tends to zeroi^ According to Fig. |31[a), one can 



see that the peaks 2, 3 and the double peaks 1 are in- 
duced by the transition |3) — >• |2) whereas the peak 4 is 
induced by the transition |2) |1). By comparing the 
energy levels and corresponding spin configurations of the 
DQD shown in Fig.[TJa) to the spectrum and DOS of the 
confined phonons in Fig. El^b) and (c) , one can also iden- 
tify the corresponding confined phonon mode and van 
Hove singularity for each peak, which are indicated by 
the arrows with the same numbers in Fig. [3][b) and (c) , 
respectively. 

It should be emphasized that not every van Hove sin- 
gularity can induce a sharp peak. This is because the 
behavior of the nominator of the coefficient T^^^^, e.g., 
the form factor, is also crucial for the peaks, which de- 
pends on the properties of the confined phonon modes 
and the electron-phonon interactions. Zeros of the form 
factors can greatly suppress the SRR. To induce a sharp 
peak, the suppression must be absent or weak at the cor- 
responding van Hove singularity. Two types of van Hove 
singularities satisfy such condition in our system. The 
first type is the van Hove singularities with q — Q for the 
radial phonon modes. The peaks (peaks 2, 3 and 4) exist 
since the corresponding form factor of the piezoelectric 
coupling tends to zero slow enough so that the suppres- 
sion effect is weak. The second type is the van Hove 
singularities with phonon wave vector q ^ 0. The corre- 
sponding peaks (the double peaks 1) are divergent due to 
the absence of zeros of the form factors at g ^ 04i Both 
the deformation potential and piezoelectric couplings can 
induce such peaks. 

These peaks can be used for electric spin relaxation 
manipulation. For example, around peak 3, by chang- 
ing the electric field within 0.01 kV/cm, the effective 
SRR can change up to 2 orders of magnitude. While 
for the electric field a little far away from peak 3, the 
effective SRRs become rather insensitive to the electric 
field. Thus it can be used as electric on-and-ofF spin 
switches. It is worth comparing our spin relaxation ma- 
nipulation scheme to the one based on GaAs DQDs in 
quantum wells where bulk phonons dominate discussed 
by Wang and Wu in Ref. There, the spin relax- 
ation manipulation relies on the fact that, for the vertical 
DQDs considered in their scheme, the electron-phonon 
scattering between the bonding and anti-bonding orbital 
states of the DQDs has a very large contribution to the 
SRR for small electric field. The SRR can be greatly 
suppressed by breaking these states by strong electric 
field, resulting in the spin relaxation manipulation. As 
the bonding/anti-bonding state plays the central role in 
the quantum-well based DQD scheme, the high energy 
anti-bonding state must be occupied or serve as the vir- 
tual state for the hyperfine-mediated spin relaxation pro- 
cess. So their scheme works with high electric field (to 
break the bonding/anti-bonding states) at high temper- 
ature (to make the high energy anti-bonding state occu- 
pied) or weak magnetic field (to ensure the hyperfine- 
mediated spin relaxation is not suppressed). In our 
scheme, the spin relaxation manipulation relies mainly 
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on the quasi-one-dimensional properties of the confined 
phonons, which does not suffer from such restrictions. It 
can also work at low temperature, strong magnetic field 
or weak electric field. Moreover, the spin relaxation ma- 
nipulation can be accomplished at different electric fields, 
which can be tuned by the nanowire radius R. Thus our 
scheme is more sensitive and flexible than the scheme in 
quantum- well-based DQDs. 

In the above discussion, facilitated with the transition 
rates, we have explained the two main features of the 
effective SRR from the equation-of-motion approach. It 
should be emphasized that in such explanation based on 
the transition rate, only the contribution of the diago- 
nal elements of the density matrix to the spin relaxation 
is considered. This is valid when the electric field is far 
away from the anti-crossing point. However, for the elec- 
tric field sufficiently close to the anti-crossing point, the 
contribution of the off-diagonal elements becomes impor- 
tant, which can have a pronounced effect on the spin 
relaxation. The spin relaxation in this case cannot be 
explained by the transition rate. This can be seen in the 
cyan dashed box in Fig.[3l^a) where overestimates the 
rate induced by the transition |3) — > |2). 

In fact, the rate at the anti-crossing point E = 
0.156 kV/cm is equal to T^2/'^^ which is just the relax- 
ation rate of the off-diagonal elements between the states 
1 3) and |2). This suggests that the off-diagonal elements 
dominate the spin relaxation at this point. This can be 
seen more clearly from the temperature dependence of 
the effective SRR in the next section. 



3. Temperature dependence 

Now we discuss the temperature dependence of the 
effective SRRs at the anti-crossing point with E = 
0.156 kV/cm and B = 1 T, shown in Fig. HJa). In the 
temperature regime we consider here, the contribution 
of the high energy levels can still be neglected, so only 
the lowest three energy levels are important. The most 
prominent feature of the effective SRRs is the existence of 
the smooth peak for f^l at T = 7 K, while is mono- 
tonically increasing. Comparing to the two rates and 
TjjT ^ , one finds that f^^ is mainly determined by the spin 
relaxation process corresponding to , which is due to 
the transition |3) — >■ |2) at the anti-crossing point. While 
fg ^ is mainly determined by which is induced by 
the transition |2) — )■ |1). The peak is mainly due to the 
transition |3) — > |2), so the peak in fj~5 is absent in f~^. 
Note that the behaviors of the time evolution of {Sz) 
from the equation-of-motion approach with the temper- 
ature before and after the peak are different, which can 
be seen in Fig. |l{b)~(e). Although both can be fitted by 
the double exponential function, {Sz) exhibits an oscilla- 
tion in the first spin relaxation process before the peak 
[Fig.|ll[b)(c)] while the oscillation vanishes after the peak 
[Fig.|ld)(e)]. 
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FIG. 4: (Color online) (a) SRRs as function of the tempera- 
ture T with E = 0.156 kV/cm and B = 1 T. Curves (dots) 
with the same color and type represent the same quantities 
in Fig. [2je). Brown solid thin curve: the relaxation rate from 
Eq. go]). Yellow dashed-dotted curve: the SRR from the 
Fermi-golden-rule approach. See text for details, (b) ~ (e) 
Time evolution of (Sz) in both the short- (lower scale) and 
long-time (upper scale) regimes with E = 0.156 kV/cm and 
B = 1 T at temperature T = 0, 3, 7 and 11 K, respectively. 
Curves with the same color and type represent the same quan- 
tities in Fig. [2] 



The above-mentioned behavior can be understood if we 
assume the off-diagonal elements of the density matrix 
between the two states |2) and |3) dominate the spin 
relaxation. As far as the two states are concerned, the 
time evolution of (Sz) due to the off-diagonal elements 
follows an oscillation decay: {Sz{t)) — A sin(ilt)e~*/^^ -I- 

B, where J7 = ^ u>fy^ — Tj"^ with representing the 

Rabi frequency and Tj"^ — F32/2 being the relaxation 
rate. The parameters A and B are constants. As the 
temperature increases, the oscillation becomes slower due 
to the increase of the transition rate F32 and vanishes 
when > uj^h, contributing another exponent decay 
to (Sz)- This explains the vanishing of the oscillation 
in (Sz) from the equation-of-motion approach at high 
temperature. 

According to the previous discussion, one can also give 
an analytical estimation for the relaxation rate induced 
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by the off-diagonal elements of the density matrix be- 
tween 1 3) and |2): 

T]^^ =T^^±e{T^^ -ooM (10) 



where VL ^ ^ T^"^ - (jO^- We plot r^^ in Fig-gfa), which 
agrees with very well. Note that the upper branch 
of [corresponding to in Eq. ((TO)) ] cannot be ob- 
tained from the equation-of-motion approach since its 
contribution to the spin relaxation is very small thus is 
difficult to be fitted numerically. 

It is noted that this peak is difficult to obtain for the 
quantum- well based QDs where bulk phonons dominate. 
This is mainly due to the fact that the transition rate 
is suppressed by the small DOS of the bulk phonons 
when the energy splitting between the two states be- 
comes small,— making it difficult to fulfill the condition 
r2"^ > ojib. Whereas the DOS of the confined phonons 
is almost constant away from the van Hove singularities, 
it is easy to have the condition T^^ > Wrb satisfied. 

It is pointed out that the large contribution of the off- 
diagonal elements makes the widely-used Fermi-golden- 
rule approachj^iiSr— which calculates the SRR by the for- 
mula Tj^^ = E,,=D,pEfif2 -^^i^^ifa' inapplicable. As an 
example, we plot the temperature dependence of the SRR 
Tq^ from the Fermi-golden-rule approach in Fig. Hl^a) 
with yellow dashed-dotted curve. It exhibits an mono- 
tonic increase without the temperature peak. More- 
over, it overestimates the SRR in the whole temperature 
regime. This suggests that for the system with strong 
spin mixing, the Fermi-golden-rule approach is inade- 
quate. To calculate the SRR properly, one has to apply 
the equation-of-motion approach. 

IV. SUMMARY 

In summary, we have investigated the effective SRRs 
induced by confined phonons in the nanowire-based 



DQDs. We first study the energies and spins of the lowest 
three states and find that both can be efficiently tuned 
by the electric field. An anti-crossing exists between the 
first and second excited states, which is due to the coac- 
tion of the electric field, the magnetic field and the SOC. 
Then we study the electric field dependence of the ef- 
fective SRR within the equation-of-motion approach and 
find that multiple sharp peaks exist in the effective SRR, 
which can be used as the electric on-and-off spin switches. 
The confined phonons are found to be crucial for these 
peaks. The quasi-one-dimensional confined phonons ex- 
hibit extremely strong DOS at the van Hove singularities, 
resulting in the multiple sharp peaks as the energy split- 
tings between states with different spins match certain 
van Hove singularities. Both the van Hove singularities 
with q 7^ and the ones with q = of radial confined 
phonon modes can induce such peaks, whilst the piezo- 
electric coupling is found to be crucial for the latter. The 
manipulation of the spin relaxation based on these mul- 
tiple sharp peaks offers on-and-off spin switches, which is 
more sensitive and flexible than the one based on DQDs 
in quantum wells. We also find that the contribution of 
the off-diagonal elements of the density matrix can dom- 
inate the spin relaxation at the anti-crossing point due 
to the large spin mixing, resulting in a smooth peak in 
the temperature dependence of the effective SRR. This 
makes the widely-used Fermi-golden-rule approach inad- 
equate in the vicinity of the anti-crossing point. These re- 
sults show a great potential of the nanowire-based DQDs 
in the spin-based quantum information processing and 
spintronic devices. 
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